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MIMO Communication System 

Bo Li, Athina R Petropulu, and Wade Trappe, 


Abstract 

Recently proposed multiple input multiple output radars based on matrix completion (MIMO-MC) 
employ sparse sampling to reduce the amount of data that need to be forwarded to the radar fusion center, 
and as such enable savings in communication power and bandwidth. This paper proposes designs that 
optimize the sharing of spectrum between a MlMO-MC radar and a communication system, so that the 
latter interferes minimally with the former. First, the communication system transmit covariance matrix 
is designed to minimize the effective interference power (EIP) to the radar receiver, while maintaining 
certain average capacity and transmit power for the communication system. Two approaches are proposed, 
namely a noncooperative and a cooperative approach, with the latter being applicable when the radar 
sampling scheme is known at the communication system. Second, a joint design of the communication 
transmit covariance matrix and the MIMO-MC radar sampling scheme is proposed, which achieves even 
further EIP reduction. 


Index Terms 

Collocated MIMO radar, matrix completion, spectrum sharing 

I. Introduction 

The operating frequency bands of communication and radar systems often overlap, causing one system 
to exert interference to the other. For example, the high UHF radar systems overlap with GSM com- 
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munication systems, and the S-band radar systems partially overlap with Long Term Evolution (LTE), 
and WiMax systems Q-Q. Spectrum sharing is a new line of work that targets at enabling radar and 
communication systems to share the spectrum efficiently by minimizing interference effects Q-Q. 
This paper investigates the problem of spectrum sharing between a MIMO communication system 


and a matrix completion (MC) based colocated MIMO radar (MIMO-MC) system |10|-112|. MIMO 


radars transmit orthogonal waveforms from their multiple transmit (TX) antennas, and their receive 
(RX) antennas forward their measurements to a fusion center for further processing. The RX antenna 
measurements could be samples of the target returns, or could be the outputs of matched filters. Based on 
the forwarded data, the fusion center populates a matrix, referred to as the “data matrix", which is then 
used by standard array processing schemes for target estimation. When the target returns are sampled 
at the Nyquist rate, and for a relatively small number of targets, the data matrix is low-rank p^ , thus, 
under certain conditions it can be reconstructed based on a small, uniformly sampled set of its entries. 
This observation is the basis of MIMO-MC radars i[I§-i|Tg, in which the RX antennas forward to the 
fusion center a small number of pseudo-randomly obtained samples of the target returns, or the result 
of matched filtering with a set of randomly selected transmit waveforms, along with information on the 
sampling scheme, with each RX antenna partially filling a row of the data matrix. Subsequently, the 
full data matrix is recovered using MC techniques. MIMO-MC radars maintain the high resolution of 
MIMO radars, while requiring significantly fewer data to be communicated to the fusion center, thus 
enabling savings in communication power and bandwidth. These savings are especially important when 
the RX antennas are on battery operated nodes, and/or the communication to the fusion center occurs in a 
wireless fashion. Compared to compressive sensing (CS) based MIMO radars, MIMO-MC radars achieve 
data reduction while avoiding the basis mismatch issues which were inherent in CS-based approaches 

m). 


In this paper, the MIMO-MC radar system is considered as the primary user of the channel, while 
the MIMO communication system is the secondary user. First, for a fixed uniformly random radar 
sub-sampling scheme, fhe communication system optimally designs its transmit covariance matrix so 
that its effective interference power (EIP) exerted to the radar RX node is minimized, while its own 
average capacity and transmit power are kept at a prescribed level. In doing so, two approaches are 
proposed, namely, a cooperative and a noncooperative approach, with the latter being applicable when 
the communication system has knowledge of the MIMO-MC radar sampling instances. It is shown that 
when the MIMO-MC radar sampling scheme is known to the communication system, the EIP can be 
greatly reduced, especially at low sub-sampling rates. Second, a joint-design of the radar sampling scheme 
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and the communication system transmit covariance matrix is proposed, targeting at minimizing the EIP 
at the radar RX node. Alternating optimization is employed to solve the optimization problem. The 
candidate sampling scheme needs to he such that the resulting data matrix can he completed. Recent 
work 0 showed that for matrix completion, the sampling locations should correspond to a binary 
matrix with large spectral gap. Since the spectral gap of a matrix is not affected by column and row 
permutations, we propose to search for the optimum sampling matrix among matrices which are row and 
column permutations of an initial sampling matrix with large spectrum gap. Even before any design is 
implemented, the MIMO-MC radar system is expected to be less susceptible to interference than a plain 
MIMO radar; this is because the interference affects only some entries of the data matrix. As it is shown 
in the paper, by appropriately designing the communication TX waveforms and/or the radar sampling 
scheme, the interference can be further reduced. 


The paper is organized as follows. Section III introduces the signal model when the MIMO-MC 
radar and communication systems coexist. The problem of a MIMO communication system sharing the 
spectrum with a MIMO-MC radar system is studied in Section IV and [V] Numerical results, discussions 
and conclusions are provided in Section Vl||VIII 


Notation: S) denotes the circularly symmetric complex Gaussian distribution with mean /r and 

covariance matrix S. | • | and Tr(-) denotes the matrix determinant and trace, respectively. The set is 
defined as L}. A/’(A) and 7^(A) denote the null and row spaces of matrix A, respectively. A*, 

and A.j respectively denote the i-th row and j-th column of matrix A. [A]jj denotes the element on the 
i-th row and j-th column of matrix A. x'^ is defined as max(0, x). 


II. Background on MIMO-MC Radars 

Consider a colocated MIMO radar system with Mt^R TX antennas and Mr^n RX antennas, targeting 
at the estimation of far-held targets. The radar operates in two phases; in the hrst phase the TX antennas 
transmit waveforms and the RX antennas receive target returns, while in the second phase, the RX 
antennas forward their measurements to a fusion center. The m-th, m G NiL antenna transmits a 
coded waveform containing L symbols ,Sm(C)} of duration Tr each. Suppose that each 

RX antenna samples the target returns with sampling interval Tr, i.e., each symbol in the waveform is 
sampled exactly once. The sampling time instances are given as {Tr, ■ • • ,LTr}. Eollowing the model 
in ug-ng, the data matrix received at the RX antennas is formulated as 

= 7pDS + Wr, (1) 
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where 7 and p respectively denote the path loss corresponding to the range hin of interest, and the 
radar transmit power; D G £,^r,RxMt,R denotes the target response matrix, which depends on the target 


reflectivity, angle of arrival and target speed (details can he found in |12|); S = [s(l), • • • ,s(L)], with 
s(/) = [si(/), • • • , SMt being the sampled waveform matrix. The transmit waveforms are typically 

orthogonal, thus it holds SS^ = I |[^. Yr = [y/?(!),... ,yi?(L)]; Wr = [wr{1), .. .,wr{L)] is the 
additive noise matrix. 

Matrix D has rank equal to the number of targets thus, it is low-rank if the number of targets is much 
smaller than Mr^R and Mt^R. Similarly, matrix DS is low-rank if the number of targets is much smaller 


than Mr^R and L. The RX antennas of the matrix completion based MIMO (MIMO-MC) radar 110|-| 121 
subsample the target returns and forward the samples, along with the corresponding sampling times to 
the fusion center, thus partially populating the data matrix. The full data matrix is then completed with 
matrix completion techniques, and target estimation can be implemented based on the completed matrix 
via standard array processing schemes 0. 


The partially filled data matrix can be mathematically expressed as follows |10|, 112| 

Clj o Yr = ftj o {'jpDS + W/j), (Scheme I) 


where o denotes Hadamard product and ftj is a matrix with “0"s or “T's, with the "l"s corresponding to 
the sampling instances. In the physical implementation, only the entries of Y^j corresponding to “l"s in 
ftj represent obtained samples. The sub-sampling rate, pj, equals \\Qi\\o/LMr^R. The above MIMO-MC 
scheme is referred to in | [T0l , fT^ as Scheme I. 

Alternatively, a random matched filter bank (RMFB) at each RX antenna generates a data matrix which 
can be expressed as GD 

fill o (YrS^) = flu o (ypD + WrS^), (Scheme II) 


where flu is a sampling matrix with binary entries and dimension Mr^R x Mt^R. The locations of “l"s 
at the m-th row are the indices of the matched filters that were used at the m-th RX antenna, i.e., 
im C The sub-sampling rate pu is defined as \\flu\\o/Mt^RM^^R. This MIMO-MC radar scheme 

is referred to in | [TT| as Scheme II. 

Early studies on matrix completion theory suggested that the low-rank matrix reconstruction from 
partial entries succeeds with high probability if the low-rank matrix satisfies the incoherence property 
and the entries are sampled uniformly at random. However, recent works | [T4| showed that, regarding 
the sampling of elements, it is sufficient that the sampling matrix has large spectral gap {i.e., large gap 
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between the largest and second largest singular values). In |10|- m that the matrix DS exhibits low 
coherence while the sampling of its elements was a result of uniformly random sampling at the RX 
antennas. 


III. System Model 

Consider a MIMO communication system which coexists with a MIMO-MC radar system as shown 
in Fig. [T] sharing the same carrier frequency. The MIMO-MC radar operates in two phases, i.e., in Phase 
1 the RX antennas obtain measurements of the target returns, and in Phase 2, the RX antennas forward 
the obtained samples to a fusion center. The communication system interferes with the radar system 
during both phases. In the following, we will address spectrum sharing during the first phase only. The 
interference during the second phase can be viewed as the interference between two communication 
systems, and addressing this problem has been covered in the literature p7| , | |T8| . 

In the following. Scheme I is used to illustrate the system model. Suppose that the two systems 


have the same symbol rate and are synchronized in sampling time (see Section VI for the mismatched 
case). We do not assume perfect carrier phase synchronization between the two systems. The data matrix 
corresponding to the radar system and the received matrix at the communication RX antennas during L 
symbol durations can be respectively expressed as 


nioYR = nio ( 7 pDS + G2XA2 + Wr), 
Yc7 = HX + pGiSAi + Wc, 


( 2 ) 

( 3 ) 


where 


• Yr, p, D, S, Wr, and Qj are defined in Secfion |II| 

. X = [x(l),... ,x(L)]; Yc = [yc(l), • • ■,yc{L)y, Wc = [wc(l),. .., wc(T)]. 

• yc(0 wc(() respecfively denofe the signal and the additive noise at the radar/communication 
RX antennas sampled at the l-th sampling time. It is assumed that wc'(/) ~ CM{0,a‘^T) and 
mvr { 1 ) ~ CAA(0,o-|l). 

• H G c^r-.cxALt.c denotes the communication channel, where and Mt^c denote respectively 


the number of RX and TX antennas of the communication system |17|; Gi G ^ denotes 

the interference channel from the radar TX antennas to the communication system RX antennas Q, 
|[^; G 2 G denotes the interference channel from the communication TX antennas 

to the radar RX antennas. It is assumed that the channels remain the same over L symbol durations. 


July 9, 2015 


DRAFT 





6 


• s(/) and x(/) respectively denote the transmit vector at the radar and the communication TX 
antennas during the /-th symbol duration. The rows of X are codewords from the code-hook of 
the communication system. 

• Ai and A 2 are diagonal matrices. The l-th diagonal entry of Ai, i.e., denotes the random 

phase offset between the MIMO-MC radar carrier and the communication receiver reference carrier 
at the l-th sampling time. The l-th diagonal entry of A 2 , i.e., denotes the random phase offset 
between the communication transmitter carrier and the MIMO-MC radar reference carrier at the l-th 
sampling time. The phase offsets result from the random phase jitters of the radar oscillator and 
the oscillator at the communication receiver Phase-Locked Loops. In the literature |[T^-||^, the 
phase jitter a{t) is modeled as a zero-mean Gaussian process. In this paper, we model {au}iLi as 
a sequence of zero-mean Gaussian random variables with variance cr^. Modern CMOS oscillators 
exhibit very low phase noise, e.g., —94 dB below the carrier power per Hz {i.e., —94dBc/Hz) at an 
offset of 27r X 1 MHz, which yields phase jitter variance cJq « 2.5 x 10“^ |22|. 

It is assumed that the MIMO channels H, Gi and G 2 are perfectly known at the communication TX 
antennas. In practice, the channel state information can be obtained through the transmission of pilot 
signals Q, |231. Based on knowledge of radar waveforms and Gi, the communication system can reject 
some interference due to the radar via subtraction. However, due to the high power of the radar Q and 
the unknown phase offset, there will still be interference in the communication received signal, i.e., 

pGiS(Ai-I)RipGiSA„, 

where Aq = diag(jaii, • • • , jaii), and the approximation is based on the fact that {aii}jL^ are small. 
The signal at the communication receiver after interference cancellation equals 


Yc = HX + pGiSA„ + Wc. 


(4) 


We observe that the residual interference is not circularly symmetric. The communication channel capacity 
is achieved by non-circularly symmetric Gaussian codewords, whose covariance and complementary 
covariance matrix are required to be designed simultaneously p4| . Here we consider the circularly 
symmetric complex Gaussian codewords x(/) ~ CJ\f{0, Rxi)> which achieve a lower bound of the channel 
capacity. The design complexity is reduced since we only need to design the transmit covariance matrix 
Rzz- 

The communication system aims at minimizing its interference to the MIMO-MC radar, while main¬ 
taining its average capacity over L symbol durations, by adapting its transmit resources in both time 
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Communication TX Communication RX 


Fig. 1. A MIMO communication system sharing spectrum with a colocated MIMO radar system 


and spatial domain. In the following two sections, the spectrum sharing problem is formulated for both 
Schemes I and II. 


IV. Spectrum Sharing with Scheme I Radars 


In this section, we design the communication transmit waveforms, and in particular their covariance 
matrix, so that we minimize the interference power at the Scheme I radar RX node, while satisfying the 
communication rate and power constraints of the communication system. The total transmit power of the 
communication TX antennas equals 


E{Tr(XX^)} = E 


Tr 




L 

^Tr(R^i), 


K \i=i J ) 1=1 

where = E{x(/)x^(/)}. 

According to Q, the total interference power (TIP) exerted at the radar RX antennas equals 

TIP = E{Tr(G2XA2A|^X^G^)} 


= EtiTr(G2R.xzG^). 


(5) 


Since the radar only forwards part of Y/j to the fusion center, only the term o (G 2 XA 2 ) represents 
effective interference to the radar system. Based on this observation, we define the effective interference 
power (EIP) at the radar RX node as 


EIP7 A E {Tr (O/ o (G2XA2) (fl/ o (G2XA2))^) } 


=E {Tr ([G2 ix(1 ) ... G2Lx(L)]A2Af [G2 ix(1) ... G2lx(L)]^) } 

=®|Tr(^EG^2ix(i)x«(/)Gi]^| 

L L 

= Y, Tr (G2iR,iG") = ^ Tr (AiG2R.;Gf), 


1=1 


1=1 


( 6 ) 
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where G 21 = ^iG 2 , with A; being a diagonal matrix whose diagonal is Ct.i, i.e., A; = diag(r2.i). We 
note that the EIP at sampling time I contains the interference corresponding to “l"s in Q.i only. It is 
equivalent to say that the effective interference channel during the /-th symbol duration is G 2 Z. 

In the coexistence model of Q and Q, both the effective interference channel G 2 Z and interference 
power at the communication receiver Rjnt/ = p‘^a‘^Gis{l)s^{l)Gi vary between sampling times. The 
communication system needs to use different covariance matrices for each symbol, i.e., in order to 
match the variation of G 2 Z and Rj^tz and minimize the effective interference to the radar system while 
maintaining the capacity. The channel can be equivalently viewed as a fast fading channel with perfect 
channel state information at both the transmitter and receiver ||25l, 1261. Similar to the definition of 


(7) 


ergodic capacity |25|, the achieved capacity is the average over L symbols, i.e., 

AC({R,z}) = ^ log 2 |I + R-;^HR,;H^| , 


where {R^z} denotes the set of all R^^z’s and R^; = Rj^tz + for all / G N^. 

In the following we will consider three spectrum sharing approaches between the communication and 
Scheme I radar, namely, a noncooperative, a cooperative and a joint design approach. In the cooperative 
and joint design approaches, the communication system knows the radar sampling scheme. The perfor¬ 
mance improvement is expected to be higher under higher level of cooperation at the cost of reduced 
security and increased coordination complexity. 


A. Noncooperative Spectrum Sharing 

In the noncooperative approach, the communication system has no knowledge of Qj. Therefore, it 
cannot obtain the expression of EIP/ of Q. In this case, the communication system will design its 
covariance matrix to minimize the TIP in Q as follows: 

(Po) ^^m ^ TIP({Ra;z}) s.t. ^^^^Tr (Ra;z) < Pt (8a) 

AC({R,z}) > C, (8b) 


where the constraint of (8a I restricts the total transmit power at the communication TX antennas to be 
no larger than P/. The constraint of ( |^ restricts the communication average capacity during L symbol 
durations to be at least C, in order to provide reliable communication and avoid service outage. {Ra;z} P 0 
imposes the positive semi-definiteness on the solution. Eet us denote by Xq the feasible set determined 
by the above three constraints. Problem (Pq) is convex. 
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The power constraints of ( [8^ and (8b I are jointly applied for all L symbol durations. The extension to 
constraints individually applied for each symbol duration is straightforward because the convexity of the 
problem is preserved. Problem (Pq) is a variant of the Problem (Pg) in | [T7| for multichannel spectrum 
sharing in cognitive radio network. 


B. Cooperative Spectrum Sharing 

In the cooperative approach, the MIMO-MC radar shares its sampling scheme Ctj with the commu¬ 
nication system. Now, the spectrum sharing problem can be formulated as 


(Pi) min EIP/({Ra.J) s.t. € Xq. (9) 

Problem (Pi) has exactly the same constraints as (Pq)- 
The Lagrangian of (Pi) can be written as 

£({R,a, Ai, Aa) =EIPK{Rxa) + As (C - AC({R,.J)) 

where Ai > 0 is the dual variable associated with the transmit power constraint, and As > 0 is the 
average capacity constraint. The dual problem of (Pi) is given as 


(Pi-D) max fl(Ai,As), 

Ai ,A2 ^0 

where ^(AijAs) is the dual function defined as 


5(Ai,As)= inf £({Ra;i}, Ai,As). 

{R2;i}^0 


The domain of the dual function, i.e., domgr, is Ai,As > 0 such that p(Ai,As) > —oo. It 
dual feasible if (Ai,As) G domp. It is interesting to note that ^(AijAs) can be obtained 
independent subproblems, each of which can be written as follows 

(Pi-sub) min Tr ((G^ A^Gs + Ail) K^i) 


is also called 
by solving L 


( 10 ) 


-Aslogs|l + R;:;/HR,zH^|. 

Before giving the solution of (Pi-sub), let us first state some observations. 

Observation 1) If there is an optimal point (and it has to be unique), the average capacity constraint is 
active at the optimal point. This means that the achieved capacity is always C and A 2 > 0. To show this, 
let us assume that the optimal point {R*;} achieves AC({R*J) > C. Then we can always shrink {R*;} 
until the average capacity reduces to C while the objective will also be reduced. Thus, we end up with 
a contradiction. 


July 9, 2015 


DRAFT 



10 


Observation 2) (G^A;G 2 + Ail) is positive definite for all I G N^. This can be shown via contradiction. 
Suppose that there exists I such that G|^A;G 2 + Ail is singular. Then it must hold that G^A;G 2 
is singular and Ai = 0. Therefore, we can always find a nonzero vecfor v lying in fhe null space 
of G 2 ^A/G 2 . At the same time, it holds that ^ 0 with very high probability, because 

H is a realization of the random channel. If we choose = avv^ and a —)• cx), the Lagrangian 
£({R 2 :z}, 0, A 2 ) will be unbounded from below, which indicates that Ai = 0 is not dual feasible. This 
means that Ai is strictly larger than 0 if G^A;G 2 is singular for any 1. The claim is proved. 

Based on the above observations, we have the following lemma. 


Lemma 1 ( |[T8|, p7|). For given feasible dual variables Ai, A 2 > 0, the optimal solution of (Pi-sub) 
is given by 


R*,(Al,A2) = ^^'^'u^5]^Uf$^'^^ 


( 11 ) 


where = G^ A;G 2 +AiI; U; is the right singular matrix o/H; = X); = diag(/3n, • ■ •, Ar) 

with j3u = (A 2 — 1/a‘ff}'^, r and au^i = 1,... ,r, respectively being the rank and the positive singular 
vales of H;. It also holds that 

r 

log 2 |I + R;1hr:,H^| = (log(A 24 ))+ ■ (12) 

i=l 

Based on Lemma [T] the solution of (Pi) can be obtained by finding fhe optimal dual variables AJ, A2. 

The cooperafive spectrum sharing problem (Pi) can be solved via the procedure outlined in Algorithm 

ffl 


Based on Lemma [T] the coexistence model can be equivalently viewed as a fast fading MIMO channel 


H;. The covariance of the waveforms transmitted on H/ is R^; = It is wcll-known that 


.1/2 


the optimum Ra^^ equals UiSjU; with power allocation obtained by the water-filling algorithm |25|. 
The achieved capacity is the average over all realization of the channel, i.e., {H;}^i- This justifies the 
definition of average capacity in (|7]). Lemma [1] shows that the communication transmitter will allocate 
more power to directions determined by the left singular vectors of H corresponding to larger eigenvalues 
and by the eigenvectors of corresponding to smaller eigenvalues. In other words, the communication 
will transmit more power in directions that convey larger signal at the communication receivers and 
smaller interferences to the MIMO-MC radars. 

The following theorem compares the minimum EIP achieved by the noncooperative and cooperative 
approaches under the same communication constraints. 
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Algorithm 1 Cooperative Spectrum Sharing (Pi) 


1: Input: Pi,C',cT^, Ai G [0,Ai] 

2: Initialization: A; = 0, A^ = Ai 
3: repeat 

4: Ai (A; + Au)/2 

5: Find the minimum A 2 > 0 such that SI=i > LC. Use the obtained A 2 to 

calculate R*;(Ai,A 2 ) according to (11). 

6: if Zti Tr(R*i(Ai, A2)) < Pt then 

7: A^ — Ai 

8 : else 

9: Xi = Ai 


10 : end if 

11: until Xu — Xi < S\, where 5^ is a predefined threshold. 
12: A); = Ai,A^ = A 2 ; 

13: Output: = R;;(A!,A^) 


Theorem 1. For any Pt and C, the ElPj achieved by the cooperative approaches in (Pi) is less or equal 
than that of the noncooperative approach via (Pq)- 

Proof: Let {R*°} and {R*)} denote the solution of (Pq) and (Pi), respectively. We know that 
{R*°} satisfies fhe consfrainfs in (Pi), which means thaf {R*°} is a feasible poinf of (Pi). The optimal 
{R*}} achieves an objective value no larger than any feasible point, including {R*^}, does. It holds that 
EIP/({R*J}) < EIP/({R*°}), which proves the claim. ■ 

There are certain scenarios in which the cooperative approach outperforms significantly the noncooper- 

1 /o 

ative one in terms of EIP. Eet us denote by the intersection of M{G 2 i) and P(Rjj H), and by (j )2 the 

1 /2 

intersection of M{G 2 ) and 'P(Rj; H). We know that (j )2 ^ fi- Consider the case where fi is nonempty 
while (j )2 is empty. This happens with high probability when Mr^R > Mtfi but pM^^n is much smaller 
than Mt^c- Problem (Pi) will guide the communication system to focus its transmission power along 
the directions in fi to satisfy both communication system constraints, while introducing zero EIP to the 
radar system. On the other hand, since 4>2 is empty. Problem (Pq) will guide the communication system 
transmit power along directions that introduce nonzero EIP. In other words, the sub-sampling procedure 
in the MIMO-MC radar may reduce the dimension of the interference channel G 2 row space. This further 
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increases the design flexibility of the communication waveforms. Therefore, it is more possible to find 
communication waveforms that satisfy the communication constraints and meanwhile introduce smaller 
EIP. 


C. Joint Communication and Radar System Design for Spectrum Sharing 

In the above described spectrum sharing strategies, the MIMO-MC radar operates with a predetermined 
pseudo random sampling scheme. However, in this section, we consider a joint design of the communi¬ 
cation system transmit covariance matrices and the MIMO-MC radar random sampling scheme, i.e., Clj. 
The candidate sampling scheme needs to ensure that the resulting data matrix can be completed. This 
means that f2/ is either a uniformly random sub-sampling matrix | [T^ , or a matrix with a large spectral 
gap 10. 

Recall that EIP/ = (^iG 2 RxzG^) . The joint design scheme is formulated as 

L 

(P 2 ) {{Rxzjjfi/} = argmin V Tr (AiG 2 Rx«Gf) 

{R.TP i=i 

s.t. {R,a G ^ 0 ? — ^iag(f^./), Cl is proper. 

The above problem is not convex. A solution can be obtained via alternating optimization. Eet ({R^^}, Cf^) 
be the variables at the n-th iteration. We alternatively solve the following two problems: 

{R-xzl = argmin V Tr (Af"^G 2 RxzGf) , (13a) 

{RxTeXo 

= argmin Tr {AiG 2 KiG^) , (13b) 

s.t. A/ = diag(ri.;), Cl is proper. 


The problem of (13a) is convex and can be solved efficiently. To avoid the intermediate variable {A;}, 


we can reformulate (13b I as 


= argmin Tr(ri^Q”) s.t. is proper, (14) 

f2 

where the l-th column of contains the diagonal entries of G 2 R";G^. Recall that the sampling matrix 
ri is proper either if it is a uniformly random sampling matrix, or it has large spectral gap. However, it 
is difficult to incorporate such conditions in the above optimization problem. 

Noticing that row and column permutation of the sampling matrix would not affect its singular values 
and thus the spectral gap, we propose to optimize the sampling scheme by permuting the rows and 
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columns of an initial sampling matrix 


ri” = argmin Tr(ri^Q”') s.t. G p(ri°), 
fi 


( 15 ) 


where denotes the set of matrices obtained hy arbitrary row and/or column permutations. The 

is generated with binary entries and [pjLMr^nl ones. Meanwhile, has large spectral gap. One of the 
matrices that exhibit large spectral gap with high probability is the uniformly random sampling matrix 
IB- Brute-force search can be used to find the optimal fi. However, the complexity is very high since 
|p(r2°)| = 0(Mr,j?!L!). By alternately optimizing w.r.t. row permutation and column permutation on 
we can solve ( [T5| ) using a sequence of linear assignment problems |28|. 

To optimize w.r.t. column permutation, we need to find the best one-to-one match between the columns 
of and the columns of We construct a cost matrix C'^ G with The 

problem turns out to be a linear assignment problem with cost matrix C^, which can be solved in 


polynomial time using the Hungarian algorithm |28|. Let denote the column-permutated sampling 
matrix after the above step. Then, we permute the rows of to optimally match the rows of Q”. 
Similarly, we construct a cost matrix G with [C'^]mi — Again, the Hungarian 

algorithm can be used to solve the row assignment problem. The above column and row permutation 
steps are alternately repeated until Tr(r2^Q’^) becomes smaller than a certain predefined threshold (5i. 

It is easy to show that the value of EIP/ decreases during the alternating iterations between ( 13a[) and 


(13bl. The proposed algorithm stops when EIP/ decreases with value smaller than a certain predefined 
threshold 62 - The proposed joint-design spectrum sharing strategy is expected to further reduce the EIP 
at the Scheme I radar RX node compared to the methods in Section IV-A and |IV-B| The complete 
joint-design spectrum share algorithm proposed in this section is summarized in Algorithm 


V. Spectrum Sharing with Scheme H MIMO-MC Radars 

When the Scheme II radar is considered, the signal model of the random matched filter can be expressed 
as follows: 

ftijo (YrS^) = ftii O (DSS^ + GsXAaS^ + WrS^). 


The effective interference power to the Scheme II radar is given by (161 on top of next page, where 


denotes the m-th row of G 2 ; Sm is composed by rows selected from S according to set as defined 
in Section |I^ Each sum term on the right hand side (RHS) of ( [T^ is the interference power at one radar 
receive antenna. To minimize the interference power with respect to the spatial spectrum we have 

the following lemma to express ([B iri terms of {R^z}- 
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Algorithm 2 Joint design based speetrum sharing between Scheme I radar and a MIMO comm, system 
1: Input: H, Gi,G2,Pt,C', 61,82 

2: Initialization: is a uniformly random sampling matrix 

3: repeat 

4: Solve problem 

5. QiPrev ^ Qn-1 

6 : loop 

7: Find the best column permutation of by solving the linear assignment problem 

with cost matrix 

8: <— Find the best row permutation of by solving the linear assignment problem with cost 

matrix 


13a I using Algorithm [f while fixing ft 


n-l 


9: if |Tr((n^)'^Q^) - < <5i then 

10 : Break 

11: end if 

12 : ^ 

13: end loop 

14: ^ 


15: n <— n + 1 

16: until |EIP” - EIP^"^| < 82 

17: Output: {R^i} = 




M, 


ElFii 4 E {Tr (nn o (GsXAsS^) {Qn o (GaXAsS^))^) } = E g^XAaS^S^Af X^g. 




m=l 


Lemma 2. For the effective interference power EIPjj, it holds that 

L 

EIPii = Y,Tr {Ai^G 2 B..iG^) , (17) 

1=1 

where = diag{ai^^,... ^); = s^(/)sm(0 Sm{l) containing entries of s{l) indexed 

by set Cm- 
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Proof: The proof can be found in Appendix ■ 

If we choose i.e., all matched filters are used and no matrix completion is considered, 

Sm equals S for all m G Then by Lemma the interference at the output of the full matched 

filter bank equals 

L 

IPpMFB — ^ a«Tr (G 2 RxzG|^) , (18) 

1=1 

where ai = s^{l)s{l). It is noted that 0 < < ai^Mm G 

In the following we discuss four levels of cooperation between the communication system and the 
Scheme II radar. 


A. Noncooperative Spectrum Sharing 

In the first case, the communication transmitter does not utilize any knowledge of the MIMO radar 


system except for the interference channel G 2 . Just as in the noncooperative case in Section IV the 
communication transmitter designs its spectrum to minimize the interference power exerted at the radar 
RX antennas, i.e., TIP, using (Pq)- 


B. Partially Cooperative Spectrum Sharing 

In the second case, the communication transmitter exploits knowledge of the a/’s, obtained by using 
shared radar waveform^ The communication transmitter designs its spectrum to minimize the interfer¬ 
ence power at the output of full matched filter banks in all the radar receivers 

(P 3 ) min IPFMFB({Ra;i}) s.t. {Rx/} G Xq. (19) 

The interference power IPfmfb has the same summation terms as in TIP but reweighed by the a/’s along 
different symbol durations. 


C. Fully Cooperative Spectrum Sharing 

In the fully cooperative case, the radar system shares the diagonal matrices A/^,Z G with the 
communication system. The spectrum sharing problem can be formulated as 

(P 4 ) min EIP//({R^./}) s.t. {Ra,-/} G Xq. (20) 

'Recall that the communication capacity in is defined based on the knowledge of S. This means that the radar waveforms 
are shared with the communication transmitter. 
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The effective interference power EIP// also has similar structure to TIP and IPfmfb in Q and 
respectively, while it is reweighed hy the diagonal matrices G We can see that the random 

matched filter hank introduces the weights A/g’s which affect the power allocation in both time and 
spatial domain. (P4) can also he solved using the dual decomposition technique used in Algorithm [T] 
The following theorem compares the effective interference power to Scheme II radar, achieved hy (Pq), 
(P3), (P4) in the above three cases. 

Theorem 2. For any Pt and C, the effective interference power to Scheme II radar achieved by (Pa) is 
not larger than those achieved by fPoj and fPsj when none or partial information is shared with the 
communication transmitter. 

Proof: Let {R*?}, {R^f} and {R*f} denote the solution of (Pq), (P3) and (P4), respectively. We 
know that both {R*°} and {R*?} satisfy the constraints in (P4), which means that {R*°} and {R*f} 
are two feasible points of (P4). Meanwhile, the optimal {R*f} achieves an objective value no larger 
than any feasible point, including {R*°} and {R*^}, does. It holds that EIP//({R*^}) < EIP//({R*°}) 
and EIP//({R*^}) < EIP//({R*f}), which prove the claim. ■ 

D. Joint Communication and Radar System Design for Spectrum Sharing 

In the above described spectrum sharing strategies, the Scheme II radar operates with a predetermined 
pseudo random sampling scheme. In this section, we consider a joint design of the communication system 
transmit covariance matrices and the MIMO-MC radar sampling scheme, i.e., Cljj. The key of applying 
the joint design scheme is to express A;^ in terms of Ctjj, which is given in the following lemma. 

Lemma 3. The effective interference power EIPjj can be equivalently expressed as 

EIPij = Tr{njMSoS)^), ( 21 ) 

where the l-th column of Q contains the diagonal entries o/G 2 Ra;zG^. 

Proof: The proof can be found in Appendix ■ 

The joint design scheme is formulated as follows 

L 

(P5) {{Rxz},$ 2 //} =argmin ^ Tr (Ai^G2RxzG|^) 
s.t. {Rxz} £ Xq, Cl G 
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As in problem (P 2 )^ a suboptimal sampling matrix ri is searched over the set of matrices obtained by 
permutating rows and/or columns of The initial Mr^R x Mt^R dimensional matrix is generated 
with [piiMt^RMr^Rl ones at uniformly random positions. This guarantees that fi® and matrices obtained 
by permutating rows and/or columns of have large spectral gap. Multiple instances of can be used 
to find a better radar sampling scheme. Similarly, the technique of alternating optimization is adopted to 
solve (P 5 ). Let be the variables at the n-th iteration. We alternatively solve the following 

two problems: 


{R”;} = argmin ^ Tr (diag (Jl" ^ (s(i) o s(/))) G 2 RX/G 

{R,„,}GXo 

r2" = argminTr(f2^Q"(SoS)^) ,s.t. n G 
o 


(22a) 

(22b) 


where the l-th column of contains the diagonal entries of G 2 R”;G^. The problem of (22al is convex 
and can be solved efficiently using Algorithm By denoting Q” = Q’^(S o S)^, subproblem (22bi can 
be formulated into exactly the same form as ( [TS] ). Analogously, ( 22 b[ ) is solved using a sequence of linear 
assignment problems p 8 |, which alternately optimize w.r.t. row permutation and column permutation on 


ri^. The corresponding cost matrices and are with entries given by [C^]mi = Q- and 

[cymi = ^LiQ?y> respectively. The complete joint-design based spectrum sharing algorithm proposed 
in this section is summarized in Algorithm 


VI. Spectrum Sharing between Mismatched Systems 


In Section III the waveform symbol duration of the radar system is assumed to match that of the 
communication system. For a typical communication channel with 40 x 10^ Hz bandwidth, the maximum 
symbol rate is 20 x 10® symbols/s. For our assumption to be valid, the radar waveform symbol duration 
need to be ^/rs, which results in a typical range resolution of 7.5 meters. 

In the following, we consider the mismatched cases. We will show that the proposed techniques 
presented in the previous sections can still be applied. Let /^ = 1 /Tr and denote the radar waveform 
symbol rate and the communication symbol rate, respectively. Also, let the length of radar waveforms 
be denoted by Lr. The number of communication symbols transmitted in the duration of Lr/J^ is 
Lc = {Lr/^ f^]. The communication average capacity and transmit power can be expressed in terms 


of {Ra;z};=i ^s in Scction IV In the following, we will only focus on the effective interference to the 
MIMO-MC radar receiver. 

If/f < fy, the interference arrived at the radar receiver will be down-sampled. Let Ii C be 
the set of indices of communication symbols that are sampled by the radar in ascending order. It holds 
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Algorithm 3 Joint design based spectrum sharing between Scheme II radar and a MIMO comm, system 
1: Input: H,Gi,G2,S,Pt,C,4,(5i,,52 

2 : Initialization: is a uniformly random sampling matrix 

3: repeat 

4: Solve problem 

5: flPrev ^ Qn-1 

6 : loop 

7: Find the best column permutation of by solving the linear assignment problem 

with cost matrix 

8 : ri'’ Find the best row permutation of by solving the linear assignment problem with cost 

matrix 


22a I using Algorithm [f while fixing ft 


n-l 


9: if |Tr((n^)'^Q^(S o S)^) - Tr((flP''^’^)^Q^(S o S)'^)| < (5i then 

10 : Break 

11: end if 

12 : ^ 

13: end loop 

14: ^ 


15: n <— n + 1 

16: until |EIP”^ -EIP77^| < 62 

17: Output: {R^i} = 


that |Xi I = Lr. Following the derivation in previous sections, we have the following interference power 
expressions: 

EIP/ = Tr (Az-GaRxzG^) , 

= Eiex, Tr (Az-jGaRx/G^) , 

where I' G is the index of I in ordered set Xi. We observe that the communication symbols indexed 
by \Xi, which are not sampled by the radar receiver, would introduce zero interference power to 
the radar system. 

If /f > f^, the interference arrived at the radar receiver will be over-sampled. One individual 
communication symbol will introduce interference to the radar system in [f^/f^\ consecutive symbol 
durations. Let Ii be the set of radar sampling time instances during the period of the f-th communication 
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k;G2Rx«G^ 


Vi^G2^xlG2 


symbol. Note that Ti is with cardinality L/s^//FJ> the collection of sets Xi,... is a partition 
of . The effective interference power for both schemes of MIMO-MC radar is respectively 

EIP" = Et° ■f'' 

where A; = Xlpex, ^i' A;g = Yhv^ii We observe that each individual communication transmit 
covariance matrix will be weighted by the sum of interference channels for radar symbol 

durations instead of one single interference channel. 

We conclude that in the above mismatched cases, the EIP expressions have the same form as those 
in the matched case except the diagonal matrices A; and To calculate the corresponding diagonal 
matrices, the communication system only needs to know the sampling time of the radar system. Therefore, 
the spectrum sharing problems in such cases can still be solved using the proposed algorithms of Sections 
IV and Further investigation will be considered as our future work. 


VII. Numerical Results 

For the simulations, we set the number of symbols to L = 32 and the noise variance to = 0.01. 
The MIMO radar system consists of colocated TX and RX antennas forming half-wavelength uniform 


linear arrays, and transmitting Gaussian orthogonal waveforms |10|. The channel H is taken to have 


independent entries, distributed as CA/’(0,1). The interference channels Gi and G 2 are generated with 
independent entries, distributed as CAA(0, cr^) and CAA(0, it|), respectively. We fix cr^ = cjI = 0.1 unless 
otherwise stated. The maximum communication transmit power is set to Pt = L (the power is normalized 
w.r.t the power of radar waveforms). The propagation path from the radar TX antennas to the radar RX 
antennas via the far-held target introduces a much more severe loss of power, 7 ^, which is set to — 30dB 
in the simulations. The transmit power of the radar antennas is fixed fo = po — WOOL / Mt^n unless 
otherwise stated, and noise in the received signal is added at SNR= 25dB. The phase jitter variance is 
taken to be = 10“^. The same uniformly random sampling scheme is adopted by the radar in both 
the noncooperative and the cooperative spectrum sharing (SS) methods. The joint-design SS method uses 


the same sampling matrix as its initial sampling matrix. The TFOCUS package \29\ is used for low-rank 

matrix completion at the radar fusion center. The communication covariance matrix is optimized according 

1 /2 

The obtained is used to generate x(/) = R^'^ randn(Mt .c, 1 )- 


to the criteria of Sections 


IV 


and 


V 


We use the EIP and MC relative recovery error as the performance metrics. The relative recovery error 
is defined as IlDS — E)S||i 7 ’/||DS||i;’ for Scheme I and IlD — D||i?/||D||i 7 ’ for Scheme II, where E)S and 


July 9 , 2015 


DRAFT 








20 


D are the completed results of DS and D, respectively. For comparison, we also implement a “selfish 
communication" scenario, where the communication system minimizes the transmit power to achieve 
certain average capacity without any concern about the interferences it exerts to the radar system. 

A. Spectrum Sharing between a Scheme I radar and a MIMO Communication System 



Fig. 2. Spectrum sharing with the Scheme I radar under different sub-sampling rates. Mt,R = 4, Mr,R = Mt,c = 8, Mr,c = 4. 



Fig. 3. Spectrum sharing with the Scheme I radar under different sub-sampling rates. Mt,R = 16, Mr,R = 32, Mt,c = Mr,c = 

4. 


1) Performance under different sub-sampling rates: There is a far-held stationary target at angle 30° 
w.r.t. the radar arrays, with target reflection coefficient equal to 0.2+0.ly. For the communication capacity 
constraint, we consider C = 12 hits/symhol. The suh-sampling rate of Scheme I radar varies from 0.2 
to 1. The following two scenarios are considered. 

In the first scenario, we use Mt^R = 4, M^^r = Mt^c = 8; -^r,c = 4. We plot the EIP results for 4 
different realizations of in Fig. |^a). For better visualization. Fig. |^b) shows the relative recovery 
errors averaged over all 4 realization of Fi*’. The cooperative spectrum sharing (SS) method (see (Pi)) 
outperforms its noncooperative counterpart (see (Pq)) in terms of both EIP and MC relative recovery 


error. As discussed in Section IV the EIP is significantly reduced by the cooperative SS method when 
p < 0.6, i.e., when pMrfi is much smaller than Mt^c- The cooperative SS method performs almost the 
same as the joint-design method in this scenario. One possible reason is that the row dimension of F2 is 
too small to generate sufficient difference in EIP among the permutations of F7. 
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In the second scenario, we choose Mt^R = 16, = 32, Mt^c = 4, Mr^c = 4. In Fig. [^a), we 

plot the EIP corresponding to 4 different realization of OP. Again, Fig. [^h) shows the relative recovery 
errors averaged over all 4 realization of The cooperative SS method outperforms the noncooperative 
SS one only marginally. This is due to the fact that both G 2 and G 2 Z are full rank. The joint-design 
method for SS in Section |IV-C| optimizes FI starting from the same sampling matrix used hy the other 
three methods. Fig. suggests that the joint-design SS method achieves smaller EIP and relative recovery 
errors than the other three methods. 

We should note that when p decreases, the null space of G 2 z expands with high prohahility, and 
the EIP of the cooperative SS method is reduced. However, if p is too small, the MC recovery at the 
fusion center fails. In the above scenarios, we would like p > 0.4 for a small relative recovery error 
in matrix completion. However, values of p > 0.6 require more samples while achieving little, or even 
no improvement on the relative recovery error. Therefore, the optimal range of p is [0.4,0.6], where 
the proposed joint-design SS method reduces the EIP by at least 20% over the “selfish communication 
method". In conclusion, the sub-sampling procedure in Scheme I radar is beneficial in terms of reducing 
the effective interference power from the communication system as well as reducing the amount of data 
to be sent to the fusion center. In addition, simulations indicate that the communication average capacity 


constraint holds with equality in both scenarios, confirming observation (1) of Section IV-B 



Fig. 4. Spectrum sharing with the Scheme I radar under different capacity constraints C. Mt,R = 4, Mr,R = Mt,c = 
8, Mr^c = 4. 


2 ) Performance under different capacity constraints: In this simulation, the constant C in the commu¬ 


nication capacity constraint of (8b) varies from 6 to 14 bits/symbol, while the sub-sampling rate p is fixed 
to 0.5. Pour different realizations of are considered. Pig. shows the results for Mt^R = 4, Mr^R = 
Mt^c = 8) Mrfi = 4. Por the “selfish communicafion" and noncooperafive SS methods, the EIP and 
relative recovery errors increase as the communication capacity increases. In contrast, the cooperative 
and joint-design SS methods achieve significantly smaller EIP and relative recovery errors under all 
values of C. This indicates that the latter two SS methods successfully allocate the communication 
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Fig. 5. Spectrum sharing with the Scheme I radar under different capacity constraints C. Mt^R = 16, M^^r = 32, Mt^c = 
Mr,C = 4. 



Fig. 6. Spectrum sharing with the Scheme I radar when multiple targets present. Mt,R = 16, Mr,R = 32, Mt^c = Mr,c = 4, 
p = 0.5 and C = 12 bits/symbol. 


transmit power in directions that result in high communication rate, hut small EIP to the Scheme I radar. 

The results for = 16, Mr^n = 32, Mt^c = = 4 are shown in Fig. Since is much 

larger than the cooperative SS method outperforms the noncooperative counterpart only marginally. 

Meanwhile, the joint-design SS method can effectively further reduce the EIP and relative recovery errors. 

3) Performance under different number of targets: In this simulation, we fix p = 0.5 and C = 12 
and evaluate the performance when multiple targets are present. The target reflection coefficients are 
designed such that the target returns have fixed power, independent of the number of targets. We observe 
that the EIPs of different methods remain constant for different number of targets. This is because the 
design of the communication waveforms is not affected by the target number. Fig. shows the results of 
the relative recovery error, which increases as the number of targets increases. All methods have large 
recovery error for large number of targets, because the retained samples are not sufficient for reliable 
matrix completion under any level of noise. The proposed joint-design SS method can work effectively 
for the Scheme I radar when a moderate number of targets are present. 
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Fig. 7. Spectrum sharing with the Scheme I radar under different levels of radar TX power. Mt^R = 16, Mr^n = 32, Mt^c = 
Mr,C = 4. 


4) Performance under different levels of radar TX power: In this simulation, we evaluate the effect of 
radar TX power p 2 , while fixing p = 0.5, C = 12 and the target number to he 1. Fig. shows the results 
of EIP and relative recovery errors for Mt^R = 16, M^^r = 32, Mtfi = = 4. Again, we see that 

the joint-design SS method performs the best, followed by the cooperative and then the noncooperative 
one. When the radar TX power increases, the EIP increases but with a much slower rate. Therefore, 
increasing radar TX power improves the relative recovery errors. 



Interference channel Gi variance ct? 



Fig. 8. Spectrum sharing with the Scheme I radar under different channel variance af for the interference channel Gi. 

Mt,R = 16, Mr,R = 32, Mt^c = Mr,c = 4. 


5 ) Performance under different interference channel strength: In this simulation, we evaluate the effect 
the interference channel Gi with different erf, while fixing p = 0.5, C = 12 and fhe targef number to 
be 1. As the communication RX gets closer to the radar TX antennas, erf gets larger. Pig. shows 
the results of EIP and relative recovery errors for Mt^R = 16, Mr^ = 32, Mtfi = Mrfi = 4. Por all 
the SS methods, when the interference channel Gi gets stronger, the communication TX increases its 
transmit power in order to satisfy the capacity constraint. Therefore, the EIP and the relative recovery 
errors increases with the variance af. We also observe that the joint-design SS method performs the best, 
followed by the cooperative and then the noncooperative one. 
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Fig. 9. Spectrum sharing with the Scheme II radar under different sub-sampling rates. Mt^R = 4, Mr,R = Mt^c = 8, Mr,c = 4. 



Fig. 10. Spectrum sharing with the Scheme II radar under different sub-sampling rates. Mt,R = 16, Mr,R = 32, Mt,c = 
Mr,C = 4. 


B. Performance of the Scheme II radar and a MIMO Communication System Spectrum Sharing 

I) Performance under different sub-sampling rates: Here we consider one far-field stationary target at 
angle 30° w.r.t. the radar arrays, with target reflection coefficient equal to For the communication 

capacity and power constraints, we consider C = 12 hits/symhol and Pt = L. Again, the suh-sampling 
rate of Scheme II radar varies from 0.2 to 1. The following two scenarios are considered. 

In the first scenario, we consider Mt^n = 4, Mr^n = Mt^c = Sj Mr^c = 4. The EIP to the Scheme II 
radar and the relative recovery errors of the matrix completion are shown in Fig. The EIPs are shown 
for four realizations of while the relative recovery errors are the average over the realizations of 
We observe that all four proposed SS methods achieve significant interference reduction compared 
to the “selfish communication method". Higher level cooperation between the MIMO-MC radar and 
communication systems achieves greater EIP reduction. The fully cooperative and joint-design SS methods 
outperform their partially cooperative, and noncooperative counterparts, which validates the statement in 
Theorem]^ However, for small p’s, the improvement achieved by the fully cooperative SS method is not 
as significant as that when the Scheme I radar is considered (see Pig. |^. This is reasonable because 
in the expression of EIP// is always full rank even for small values of p. Decreasing p will not reduce 
the rank of effective reference channel Therefore, the communication system cannot find a 

direction that would introduce zero EIP to the radar. 
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In the second scenario, we consider Mi^r = 16, M^^r = 32, Mt^c = Mr^c = 4. Fig. 10 shows the EIP 
and the relative recovery errors of the matrix completion. Again, the EIPs are shown for four realizations 
of f2°, while the relative recovery errors are the average over the realizations of The joint-design 
SS method achieves much smaller EIP and relative recovery errors than the other four methods. This 
validates the effectiveness of the proposed joint-design SS method for the Scheme II radar. We conclude 
that the MC approach henefits the Scheme II radar hy reducing not only the data to he forwarded to the 
fusion center hut also the effective interference from the communication system when spectrum sharing 
is considered. 



—Selfish Communication 
— Noncooperative SS 
-Partially cooperative SS 
-Fully cooperative SS V’’’ 

—Joint-Design SS _ 


Selfish Communication 
Noncooperative SS 
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Fully cooperative SS 
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Capacity Constraint C 



8 10 
Capacity Constraint C 


Fig. 11. Spectrum sharing with the Scheme II radar under different capacity constraints C. Mt,R = i, Mr,R = Mt,c = 
8, Mr,c = 4. 
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Fig. 12. Spectrum sharing with the Scheme II radar under different capacity constraints C. Mt,R = 16, Mr,R = 32, Mt,c = 
Mr,C = 4. 


2 ) Performance under different capacity constraints: In this simulation, the Scheme II radar has fixed 


suh-sampling rate p = 0.5, while the communication capacity constant C in (8hi varies from 6 to 14 
hits/symhol. Similarly, two scenarios are considered. The results for the scenario of Mt^R = 4, Mj-^r = 


Mt^c = 8) = 4 are plotted in Pig. 11 and the scenario of Mt^R = 16, Mr^R = 32, Mt^c = = 4 


in Pig. 12 The “selfish communication" mefhod is inferior fo all proposed SS mefhods. One can also 
observe fhat fhe joinf-design SS mefhod always achieves considerably smaller EIP and relafive recovery 
errors fhan fhe partially and noncooperafive SS mefhods. 
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Fig. 13. Spectrum sharing with the Scheme II radar when multiple targets present. Mt,R = 16, Mr,R = 32, Mt,c = Mr,c = 4, 
p = 0.5 and C = 12 bits/symbol. 


3) Performance under different number of targets: In this simulation, we fix p = 0.5 and C = 12 and 
evaluate the performance when multiple targets are present. The target reflection coefficients are designed 
such that the target returns have fixed power, independent the number of targets. Again, we observe that 
the EIPs of different methods remain constant for different number of targets. The results of the relative 


recovery error are shown in Fig. 13 The proposed joint-design SS method can work effectively for the 


Scheme II radar when a moderate number of targets are present and sufficient samples are used for matrix 
completion. 


Selfish Communication 
Noncooperative SS 
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Fig. 14. Spectrum sharing with the Scheme II radar under different levels of radar TX power. Mt,R = 16, Mr,R = 32, Mt,c 
Mr,C = 4. 


4) Performance under different levels of radar TX power: In this simulation, we evaluate the effect 


of radar TX power p 2 , while fixing p = 0.5, C = 12 and the target number to be 1. Fig. 14 shows the 


results of EIP and relative recovery errors for Mt^R = 16, Mr^R = 32, Mt^c = = 4. We can see 

that the joint-design SS method greatly outperforms the other three methods. When the radar TX power 
increases, the performance gap between the joint-design SS method and the other three methods becomes 
larger. 
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Fig. 15. Spectrum sharing with the Scheme II radar under different channel variance af for the interference channel Gi. 

Mt,R = 16, Mr,R = 32, Mt,c = Mr,c = 4. 


5 ) Performance under different interference channel strength: In this simulation, we evaluate the effect 
the interference channel Gi with different af, while fixing p = 0.5, (7 = 12 and the target number to 


be 1. Fig. 15 shows the results of EIP and relative recovery errors for Mt^R = 16, Mr^R = 32, Mt^c = 
Mr,c = 4. For all the SS methods, when the interference channel Gi gets stronger, the communication 
TX increases its transmit power in order to satisfy the capacity constraint. Therefore, the EIP and the 
relative recovery errors increases with the variance af. We also observe that the joint-design SS method 
greatly outperforms the other three methods. 


VIII. Conclusions 

This paper has considered spectrum sharing (SS) between a MIMO communication system and a 
MIMO-MC radar system using two different schemes. Scheme I and Scheme II. In order to reduce the 
effective interference power (EIP) at radar RX antennas, we have first proposed two communication 
transmit covariance matrix design strategies, namely, a noncooperative and a cooperative SS method, for 
both Schemes I and II. Our theoretical results guarantee that the cooperative approach can effectively 
reduce the EIP to a larger extent as compared to the noncooperative approach. Second, we have proposed 
a joint design of the communication transmit covariance matrix and the radar sampling scheme to further 
reduce the EIP. The EIP reduction and the matrix completion recovery errors have been evaluated under 
various system parameters. We have shown that both Scheme I and II radars enjoy reduced interference 
by the communication system when the proposed SS methods are considered. In particular for Scheme 
I, the sparse sampling at the radar RX antennas can reduce the rank of the interference channel. Our 
simulations have confirmed that significant EIP reduction is achieved by the cooperative approach; this 
is because in that approach, the communication power is allocated to directions in the null space of 
the effective interference channel. When the number of radar RX antennas is much larger than that 
of the communication TX antennas, the cooperative approach outperforms the noncooperative one only 
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marginally. Our simulations have suggested that for both Schemes I and II, the joint-design SS method 
can achieve much smaller EIP and relative recovery errors than other methods when the number of radar 
TX and RX antennas is moderately large. 


Appendix A 
Proof of Lemma[2] 

Let us first look at the interference power at the m-th radar receive antenna, m G 

= ^ {g^XAsS^S^A^X^g^} 

= E {Tr (g^XAsS^S^A^X^g™)} 

= E{Tr(s|^A^X^ g^XAa) } ^23) 

=A„ 

= Tr (A^E {Al^X^B^XAs}) 

^Tr(A™C™), 

The entry on the fe-th row and l-th column of Cm equals 

= E |eJ(“^'-“2'‘)Tr (x^(/c)Bmx(0)} 

= E |eJ(“='-“='')Tr (B^x(0x'^(fe))} 

= e^'("=‘-“^'=)Tr (B^E {x(f)x^(A;)}) 

= f^fciTr (BmRxi) ; 

where Sm is the Kronecker delta function with value 1 if A; = f, and value 0 otherwise. Therefore, matrix 
Cm is diagonal with Tr(BmRxi) as its f-th entry. The interference power at the m-th radar receive 
antenna can be expressed as 

13m = Tr (A^diag (Tr (B^R^i), ■ • •, Tr (BmR 2 ,.L))) 

= ai5„Tr(BmRa:z) 

(24) 

= Xl/=i (gmRxigm) 

= ai^^Sm'^xlgm 
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where denotes the l-th diagonal entry of Am as defined in (jlTj). Substituting f3m in (241 into ( [T^ , 
we obtain the expression of the effective interference power to the Scheme II radar as follows 






EIP/7 ='^ Pm =Y^ 

m=l m=l 1=1 

L Mr-m 

1=1 m=l 
L 

= J]Tr(Az^G2R:./G|^), 


Z=1 


which completes the proof. 


Appendix B 
Proof of Lemma[3] 

The m-th diagonal element of can be expressed as 

= s^{l)A^^s{l) = Tr (s^(f)A^^s(f)) 


(25) 


= Tr (A^^s(f)s^(f)) = Qm- (s(0 o s(/)), 

where Ag^ = diag(nm ) and Qm- denotes the m-th row of flu. Thus we have A;^ = diag {flu (s(f) o s(Z))). 
Substituting this into ( [TT] ) gives 

L 

EIP// = ^ Tr (diag (f2// (s(Z) o s(Z))) G 2 RX 1 G 2 ) 

1=1 

= Tr {[n 7 /(s(l) o s(l)),..., nn{s{L) o s{L))f q} 

= Tr {[12/7(8 o S)]^Q} = Tr o Sf) . 

Lemma is proved. 
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